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1. Introduction 

Detached eclipsing binaries provide precise fundamental stel- 
lar parameters like mass and radius and are thus the prerequi- 
site for the validation of stellar evolutionary models. The em- 
pirical constraints from over four dozen systems have shown 
that for main sequence stars between 1 and lOM^ the agree- 
ment is acceptable, i.e. better than 2% dAndersenll99ll [l998). 
while at the lower main sequence the situation is far less satis- 
fying. Up to now, only three eclipsing systems with M-type pri- 
maries are known, despite the fact that low mass main sequence 
stars dominate the stellar population by number. The first such 
system to be discovered was YYGem (Jov & Sanford 1926; 
Ivan GenJl926T) . followed by CM Dr a dEggen & Sandagel 19671 
and references therein) and CU Cnc (Delfosse et al. 1999); first 
mass determinatio ns came from iLeung & Schneider! Q1978), 
Lacvl dl977i and belfosse et all fl999l) . respectively. While 
Metcalfe et alJ d 19961) find the slope of the mass-radius rela- 
tion derived from the M dwarf binary system C M Dra in agree- 
ment with model predictions, iDelfosse et ail feOOOl) reported 
on a disagreement between empirical and theoretical mass- 
luminosity relations of 10-2 0% in the V band, a nd recent pre- 
cise analyses of YY Gem tTorr es & Ribasl 120021) and CU Cnc 
(Ribas l2003h also revealed an underestimation (10-20%) of 
the radii of low mass stars from current evolutionary mod- 
els. Additional constraints for the empirical mass-radius rela- 
tion are provided by the first interferometr ic measurements o f 
radii from lower main sequence stars ( Segrans an et alJ l2003). 
These results agree well with model predictions at the present 
accuracy level, with a possible discrepancy for stars with 0.5- 
0.8Mg. Such observations do not provide an independent 
measurement of the stellar mass, however, so that eclipsing 
systems still are the primary source for a model-independent 
determination of fundamental parameters. 

Future improvements of the theoretical mass-radius re- 
lation for the lower main sequence would strongly bene- 
fit from a larger empirical database through an increased 
sample of eclipsing binaries. Recently, 137 eclipsing low- 
luminosity candidates were announced by the OGLE (Optical 
Gravitational Lensing Experiment) consortium ( Ud alski et alJ 
2002albl E003l) . of w hich several of the s econdaries turned out 
to be M-type stars dPreizler et al.ll2002l) . In this paper we re- 
port the discovery of another interesting eclipsing binary sys- 
tem, 2MASSJ05 16288+260738, whose components appear to 
bracket the M-star range, with the potential of extending the 
empirical mass-radius relation into the sub-stellar range. 

The new eclipsing system has been discovered in observa- 
tional data taken during a coordinated photometric monitoring 
campaign in December 2001. This dataset has been obtained 
to monitor the light variations of the DB variable white dwarf 
KUV 05134+2605 ddraueret al.ll989i Handler et al. in prep.). 
It consists of many individual light curves taken by either pho- 
tomultiplier (PMT) or CCD instruments; the newly discovered 
object is included in 48 individual time series of images ob- 
tained with CCD cameras. While analysing field stars for pho- 
tometric stability to check whether they could be used as refer- 
ences, an object located a little less than one arcminute south 
of the DB was found to show the signature of an eclipse in 



Fig.l. Finding charts for 2MASSJ05 16288+260738 (DSS-2 
red: left, DSS-2 blue: right). The side length is 4'x4' for each 
image; north is up and east is to the right. 

the Calar Alto 2.2 m data set of 2001 Dec 07 (see Table GJ. 
Subsequent searches in the other data sets revealed that eight 
more eclipses had partly or fully been observed. A year later, 
5 additional data sets were obtained, two of which covered the 
eclipse. The full time-resolved photometric data are compiled 
in Section |3 Archive searches contributed an identification of 
the object as well as additional colour information (Sections 
|2] and |4j. Two months after the initial observations, an opti- 
cal spectrum could be obtained, and in the following observing 
season, an infrared spectrum was taken (see Section^. 

In the following, we compile the information that is cur- 
rently available on the object, report our results from the light 
curve solution and the spectral analysis, and propose a possible 
configuration for this system. 

2. Positional information 

A search with SIMBAD yielded no catalogued object at or 
near the coordinates of the eclipsing object, but loading the 
Incremental Release Extended Source Catalog of the Two 
Micron All Sky Survey (2MASS) into ALADIN resulted 
in a match. We could clearly identify our object with the 
point source 2MASS J05 16288+260738, and later with a point 
source in the USNO-B catalogue (cf. Sect. l4.1> . We use the 
2MASS catalog entry to give improved coordinates: 

RA = 05 h 16 m 28181, 5 = +26°07'38'.'8 (J2000). 
For a clear identification, the object is marked with horizontal 
bars in the finding charts given in Figure[T] 

3. The light curve 

3. 1. Time-resolved photometric data 

All photometric data sets obtained with CCD cameras in the 
December 2001 KUV 05134+2605 campaign were compiled. 
Additional observations obtained in November 2002 were 
added later. A list of all data sets used is given in Tabled with 
a complementary key to the observing sites involved in Table|2] 
All data were bias and flatfield corrected according to standard 
routines. Aperture photometry w as performed on a ll of these 
frames using the TRIPP package dSchuh et aljE 003). Two ref- 
erence stars that are available on all frames (shown to be stable 
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Table 1. Photometric observations. 



Site 


start time [UT] 




length[h] 


frames 


epoch 


minimum [HJD] 


O-C 


comments 


SARA 


2001 Dec 06 


05:01 


7.1 


857 








used for profile 


SARA 


2001 Dec 07 


03:14 


5.5 


664 








used for profile 


SARA 


2001 Dec 07 


10:14 


1.9 


233 










CAHA 


2001 Dec 07 


20:09 


8.8 


1032 





2452251.5164 


-0.0009 


used for profile 


Piszkesteto 


2001 Dec 07 


19:23 


7.3 


380 





2452251.5164 


-0.0009 




BO AO 


2001 Dec 08 


11:20 


9.0 


1029 










Piszkesteto 


2001 Dec 08 


18:14 


4.7 


200 








used for profile 


WISE 


2001 Dec 09 


01:12 


2.4 


690 








used for profile 


WISE 


2001 Dec 09 


18:02 


1.8 


511 










Piszkesteto 


2001 Dec 09 


18:26 


7.5 


370 








used for profile 


WISE 


2001 Dec 09 


20:46 


2.6 


768 








used for profile 


WISE 


2001 Dec 09 


23:43 


2.0 


564 








used for profile 


WISE 


2001 Dec 10 


01:47 


1.8 


505 








used for profile 


BO AO 


2001 Dec 11 


16:51 


1.6 


181 








used for profile 


BO AO 


2001 Dec 11 


19:17 


0.9 


101 








used for profile 


SAAO 


2001 Dec 11 


22:22 


3.0 


1080 


3 


2452255.4007 


+0.0016 




GAO 


2001 Dec 12 


11:33 


1.4 


308 










SAAO 


2001 Dec 12 


21:41 


1.7 


627 










GAO 


2001 Dec 13 


13:47 


5.3 


1000 










SAAO 


2001 Dec 13 


19:57 


5.4 


1821 










GAO 


2001 Dec 14 


12:37 


6.7 


1301 


5 








BO AO 


2001 Dec 14 


16:35 


2.5 


257 








used for profile 


Piszkesteto 


2001 Dec 14 


18:08 


6.7 


328 








used for profile 


SAAO 


2001 Dec 14 


20:11 


5.1 


1821 










SAAO 


2001 Dec 15 


22:01 


1.9 


643 










GAO 


2001 Dec 16 


16:25 


3.2 


700 








used for profile 


SAAO 


2001 Dec 16 


19:58 


5.3 


1901 


7 








GAO 


2001 Dec 17 


10:51 


0.9 


203 










GAO 


2001 Dec 17 


13:19 


0.6 


125 










GAO 


2001 Dec 17 


15:50 


3.2 


700 










SAAO 


2001 Dec 17 


20:02 


5.1 


1827 










GAO 


2001 Dec 18 


10:22 


8.6 


1750 










SAAO 


2001 Dec 18 


20:03 


4.8 


1710 










GAO 


2001 Dec 19 


12:22 


6.7 


753 


9 


2452263.1615 


-0.0013 


used for profile 


SAAO 


2001 Dec 19 


20:04 


4.8 


1738 










SARA 


2001 Dec 19 


04:54 


7.2 


817 








used for profile 


GAO 


2001 Dec 20 


10:54 


6.2 


650 








used for profile 


SAAO 


2001 Dec 20 


19:58 


4.9 


1752 


10 


2452264.4599 


+0.0031 




SARA 


2001 Dec 20 


02:37 


9.7 


1107 








used for profile 


SAAO 


2001 Dec 21 


19:53 


5.0 


1804 










GAO 


2001 Dec 22 


12:50 


6.0 


625 








used for profile 


SAAO 


2001 Dec 22 


21:22 


3.0 


1074 










SARA 


2001 Dec 22 


04:31 


3.3 


404 


11 


2452265.7487 


-0.0019 




SARA 


2001 Dec 23 


05:12 


0.8 


100 










SARA 


2001 Dec 23 


10:53 


1.0 


125 


12 








GAO 


2001 Dec 23 


12:37 


5.8 


700 


12 


2452267.0450 


+0.0004 


used for profile 


SAAO 


2001 Dec 23 


19:54 


4.9 


1663 










GAO 


2001 Dec 24 


12:37 


6.0 


720 








used for profile 


CAHA II 


2002 Oct 31 


23:11 


3.6 


50 








used for profile 


CAHA II 


2002 Nov 04 


22:43 


1.8 


27 








used for profile 


CAHA II 


2002 Nov 07 


02:14 


2.4 


17 


259 


2452586.6506 


+0.0010 




CAHA II 


2002 Nov 09 


23:02 


6.9 


99 








Johnson I filter data 


CAHA II 


2002 Nov 11 


22:59 


5.2 


75 


262 


2452590.5304 


-0.0010 


Johnson I filter data 



Table 2. Key to observatory sites. 



Site Telescope Observers 



CAHA 


Calar Alto Observatory, Centro Astronomico Hispano Aleman, Almerfa, Spain 


2.2 m 


SD, SLS 


CAHA II 


Calar Alto Observatory, Centro Astronomico Hispano Aleman, Almerfa, Spain 


1.2m 


TN, EG 


BO AO 


Bohyunsan Optical Astronomy Observatory, Korea 


1.8m 


SLK, BGP 


SAAO 


South African Astronomical Observatory, Sutherland, South Africa 


1.0m 


GH, TM 


SARA 


Kitt Peak National Observatory, Tucson, Arizona, United States of America 


0.9 m 


MW 


GAO 


Gunma Astronomical Observatory, Japan 


1.5m 


OH, KK, HT, EK 


WISE 


The Florence and George Wise Observatory, Tel-Aviv University, Israel 


1.0m 


EL, PI, YL 


Piszkesteto 


Piszkesteto, the mountain station of Konkoly Observatory, Matra Mountains, Hungary 


1.0m 


MP, BS, GV, DZ 
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during the whole 2001 campaign) were chosen and used con- 
sistently for all data sets to produce relative light curves. The 
total light curve was then scaled by a unique factor to produce a 
light curve with a mean relative intensity of unity for the white 
light contributions outside eclipse. Finally, all times were con- 
verted from Julian date (JD) to heliocentrically corrected Julian 
date (HJD). The result is shown in Fig.|2] 

The light curve shows a clear periodicity of 1.29 days. All 
observed eclipses are similar to each other, have a duration 
of about 0.10 days and exhibit a decrease in flux of 15% (or 
0.17 mag) at the deepest point. There is no indication of a sec- 
ondary eclipse in any of the eight 2001 data sets that partly or 
fully cover the phase where such an event would be expected. 
Furthermore, the three 2002 data sets covering that phase put a 
clear upper limit to the depth of any secondary eclipse: 0.49 % 
(or 5.4mmag) in white light and 0.70% (or 7.6mmag) in the 
Johnson I filter. 

3.2. Ephemeris 

Primary minima times were determined by fitting parabolas to 
the eclipses. The results for the epochs (two independent data 
sets), 3, 9, 10, 11, 12 (concatenated from two different non- 
overlapping data sets), 259 and 262 are given in Tabled no fits 
could be obtained for epochs 5 and 7 since only parts of either 
ingress or egress had been observed there. A linear regression 
for the measured minima times then gives the linear elements 
and their 1 a errors for the primary minima as 

HJD = 2452251 d 5173 + l d 29395 • E. 

±16 ±25 

This ephemeris was used to generate a folded profile from the 
data taken in 2001. 

The folded profile has also been carefully inspected to ver- 
ify that no secondary eclipse is apparent in the data. The profile 
remains at the same relative flux level outside of the primary 
eclipse with no significant indication of ellipsoidal light vari- 
ations or reflection effects. It was then used to obtain a light 
curve solution as discussed in Sec.|6] 

4. Colour and spectral information 

4.1. Colours 

To derive a visual magnitude for the analysed object, we have 
obtained further photometry of the field. The observations 
were made at the South African Astronomical Observatory's 
30" telescope using the UCT CCD Photometer. A series of 
10 images was taken in a Johnson V filter on March 10 2003 
starting at 18:19:03 UT (exposure time 100 s). The measure- 
ments were thus made out of eclipse. We compare our results 
for KUV 05134+2605 and 2MASS J05 16288+260738 to esti- 
mate that the mean V magnitude difference (KUV — 2MASS) 
is 1 .44 m ±0.02. Due to the rapid nature of the variations in the 
DB variable, it is not a problem to der ive a good mean m agni- 
tude for that object. Using V=16.70 (Wegner et"aTl l 19901) 1 for 
KUV 05134+2605, we obtain V 2M ass = 18.1 m ±0.1. 



Photographic B, R and I mag nitudes are published in the 
USNO-B catalogue i tMonet et alJl2003l) . and infrared J, H and 
K magnitudes are available from the 2MASS catalog. They are 
compiled in TablefJ The USNO-B and 2MASS magnitudes m 
have been converted to F\ using the relation 



Table 3. Johnson V, USNO-B B, R, I and 2MASS J, H, K mag- 
nitudes. 







B 


V 


R I 


j 


H K 


maj 


;nitude 


19.47 


18.1 


16.8 15.84 


14.247 


13.346 13.115 


error 


±0.3 


±0.1 


±0.3 ±0.3 


±0.040 


±0.039 ±0.039 


Ao 


H 


0.43 


0.55 


0.70 0.90 


1.25 


1.65 2.17 




[Jy] 


4440 


3810 


2880 2240 


1593 


1089 713 




[ * ] 


1.17 


2.17 


3.36 3.83 


6.12 


5.51 2.58 


* 


[erg s" 


1 cm _2 A" 


- 1 ■ 10 


-16] 








ergs - 


W 2 A 




i^[Jy]-i0" 


0.4m _ g 


1Q- 13 /A [/H 2 



The photometric zero points F Va and central wavelengths 
Ao used for the conversion are tabulated in TablefJ] along with 
the results for F\ . For the 2MASS fi lters, these quantities were 
obtained from ISquires et all ll2002t) . w hile for the whi t e light 
we used the values pu blished for V by Campi ns et alJ £l985) 
and lRieke etail J 19851) for the Johnson UB VRI+ system. 

The Sloan Digital Sky Survey (SP SS) does not cover the 
field in its EDR (Early Data Release, Stought on et alJl2002l) 
so that no further photometric information is available. Since 
the object, according to its infrared colours, is very red, we 
also checked the VLA FIRST survey at 20 cm, which cur- 
rently does not cover this field eith er, and the NRAO/VLA 
Sky Survey (NVSS) at 1.4 GHz JCondon et alJll998t) . which 
covers the field but does not show a radio source in the vicin- 
ity. For completeness, we finally note that neither the ROSAT 
Bright So urce Catalogue as compiled from the WFC All Sky 
Survey JPounds et alJ Il993h nor the ROSAT XUV Pointed 
Phase Source Catalogu e as compiled from WFC observations 
during pointed phase JKrevsing et aljfl995l) list sources at or 
near the object's position. 

4.2. Optical spectra 

Two medium resolution spectra of 2MASSJ05 16288+260738 
were obtained in February 2002 at the Calar Alto 3.5 m tele- 
scope with the double beam spectrograph TWIN (see Tabled 
first part). Gratings #5 and #6 were used for the blue and red 
arm, respectively, with the dichroic set at 550 nm. Together 
with slit widths of l'/2 and l'/5 for the first and second ex- 
posure, this resulted in spectral resolutions of 0.94 and 1 .04 A. 
Both spectra turned out later on to have been taken well out- 
side any eclipse, but only the spectrum taken on February 25 
reaches an exposure level acceptable for further analysis: the 
signal-to-noise level per pixel for the first spectrum is only 3, 
but reaches 8 for the second one. The frame and corresponding 
wavelength calibration frame were bias and flatfield corrected. 
Then the spectrum was extracted, sky corrected, subjected to 
a cosmic ray filtering, corrected for the illumination function, 
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Fig. 2. Overview of the photometric observations; the flux is given in fractional intensity units. Time is in days, with data binned 
into units of 30 s. Time increases from left to right and from bottom to top. The epochs (labelled #n) are displayed continuously 
up to the end of the 2001 campaign, while for the 2002 observations only those epochs were included in the plot for which data 
points exist. The primary eclipse is displayed at multiples of 0.97 d to place it at a phase of 0.75 in this plot, allowing convenient 
viewing of both phase 0.25 where the secondary minima would be located as well as of the primary eclipses (both marked by 
horizontal dashed lines). The scatter in the individual light curves contributed by different sites varies according to aperture, 
actual exposure time and weather conditions. 
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Table 4. Spectroscopic observations 



object 


instrument 


A range [A] 




start time (UT) 




exptime 


HJD 


2MASS J05 16288+260738 
G191-B2B 

2MASS J05 16288+260738 


TWIN 
TWIN 
TWIN 


3900 
3900 
3900 


-5000, 6000 
-5000, 6000 
-5000, 6000 


-7090 
-7090 
-7090 


2002 Feb 23 
2002 Feb 25 
2002 Feb 25 


21:34 
18:40 
20:49 


1800 s 
300 s 
1800 s 


2452329.411 
2452331.281 
2452331.380 


2MASS J05 16288+260738 
GD71 


OMEGA-Cass 
OMEGA-Cass 


HK: 
HK: 


14000- 
14000- 


25000 
25000 


2002 Oct 27 
2002 Oct 28 


23:31 
02:21 


24 x 120 s 
10 x 120 s 


2452575.484 
2452575.602 


2MASS J05 16288+260738 
GD71 


OMEGA-Cass 
OMEGA-Cass 


JH: 
JH: 


10000- 
10000- 


18000 
18000 


2003 Feb 03 
2003 Feb 03 


21:49 
23:12 


24 x 120 s 
24 x 120 s 


2452674.412 
2452674.470 



and finally wavelength calibrated. Flux calibration was done by 
first applying the same procedure to an exposure of the stan- 
dard star G191-B2B taken in the same night, then using tabu- 
lated flux values to do the absolute calibration. The two result- 
ing optical spectra in the wavelength ranges of 3900 -5000 A 
and 6000-7090A of 2MASSJ05 16288+260738 are displayed 
in Fig.|5](rescaling as described in Sect.|5}- At the given S/N 
ratio, no lines or features, in particular no TiO bands, can be 
discerned. 

4.3. Infrared spectra 

By the start of the following observing season for 
2MASSJ05 16288+260738 in autumn 2002, the compila- 
tion and reduction of the light curve was not only complete 
enough to allow the prediction of eclipse times, but also to 
attempt a first light curve solution, based on an estimate of the 
spectral class obtained from the slope of the optical spectrum. 
This confirmed the suspicion that the system might be made 
up from two low mass stars or a low mass star and a substellar 
object, and therefore justified taking infrared spectra during 
Director's discretionary time at Calar Alto Observatory. An H 
and K band spectrum was observed in October 2002, and a J 
and H band spectrum in February 2003. 

Each time, a set of 24 spectra was obtained using the 
OMEGA-Cass instrument mounted on the 3.5 m telescope (see 
Table 0), well off both primary and secondary eclipse. Since 
the background is high for infrared observations, the set of 24 
spectra was obtained in such a way that alternating exposures 
contain the source on two different locations on the chip. After 
bias and flatfield correction of the individual exposures, this al- 
lows to determine a mean background at the (dispersed) loca- 
tion of the source for both types of images from the respective 
subset of the complementary frames. These two measures of 
the mean background can then be used to subtract the appro- 
priate background from all frames of the two subsets. To do 
this, the shift for each background row was first determined by 
cross-correlating it along the dispersion direction with the cor- 
responding image row, and the overall run of the shift along the 
chip obtained by fitting these row-by-row measurements with 
a low-order polynomial. To achieve the best possible subtrac- 
tion in the vicinity of the source location on the chip, this fit 
with sub-pixel accuracy was then used to shift the background 
onto the image before subtracting it. For the wavelength dis- 
persions at the two source locations, no shifts could be detected 



during the course of the exposure series. Therefore, next these 
bias subtracted, flatfield corrected and background subtracted 
frames of each of the two sets were added to yield two summed 
images. The spectra were extracted from these two images us- 
ing standard procedures for extraction, sky correction, cosmic 
filtering, illumination correction, and wavelength calibration. 
The same procedure was used for the set of 10 and 24 expo- 
sures of the standard star GD 71. For each observation, the two 
spectra for the standard star were then combined and compared 
to tabulated flux values to obtain the factor for absolute flux cal- 
ibration, which was then applied to both the two object spectra 
and the two individual standard star spectra. 

Comparing the results for the individual spectra for both 
stars leads to the conclusion that the error bars in the resulting 
combined spectra must be considered to be of the same mag- 
nitude as any "features" that one might be tempted to spot.The 
same conclusion results from a comparison of the H band parts 
of the spectra from the two different observing dates, where 
most "features" are not reproduced. Furthermore, a flux differ- 
ence by a factor of about 1.5 between those two independent 
observations gives an estimate for the errors in the flux cali- 
bration. The rescaled infrared spectra for 14000-25000 A and 
10000- 18000 A are displayed in Fig.|3] 

5. Spectral analysis 

The flux calibrated optical and infrared spectra as well as the 
broad band filter measurements converted to flux values (dia- 
mond symbols) are all displayed together in Fig. [3] To obtain a 
consistent image of the spectral energy distribution, the mag- 
nitude measurements were used to rescale the spectra where 
necessary. A unique correction factor was applied to both opti- 
cal spectra simultaneously, and a correction factor was applied 
to each of the two independent infrared spectra (HK and JH). 

The uncertainty in the optical spectrum results since 
both object and flux standard star were observed under non- 
photometric conditions. The same argument applies to the in- 
frared spectral observations, where observations from different 
nights, although both nominally flux calibrated, result in differ- 
ent flux levels for the overlapping H band. A consistent adjust- 
ment therefore seems justified. Residual errors may result from 
the transformation of magnitudes to Fa- 

In the following, it will be assumed that the observed spec- 
tral energy distribution consists of light from the primary only; 
furthermore, for reasons detailed in Sect. 17.21 the primary will 
be presumed to be a late main sequence star. 
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Since 2MASS J05 16288+260738 is located close to the 
galactic plane, the effect of interstellar reddening is not neg- 
ligible even for low-luminosity and close-by objects. 

For the initial analysis of the observed data we use a grid of 
model atm ospheres and synt hetic spectra that is based on the 
models of lAllard et alJ J200ll) . We have extended the model 
grid to effective temperatures of 10000K for gravities from 
5.5 < \ogg < —0.5 using spherical symmetry. The mix- 
ing length was set to twice the pressure scale height, this 
choice of the mixin g length was calibrated on early M dwarfs 
jLudwig et alfcOOa . 

Synthetic spectra generated from the models were com- 
pared to the observed spectra using an IDL program. This step 
was restricted to the infrared spectra. First, the resolution of the 
synthetic spectra was degraded to that of each observed spec- 
trum by convolution with a Gaussian of the appropriate width, 
and the spectra were normalized to unit area for scaling. Next, 
for each observed spectrum the program calculated a quality 
function q, similar to a \ 2 value, for the comparison with all 
synthetic spectra in the grid. The quality function is calculated 
by first scaling the model spectrum to the observed fluxes and 
then by mapping the synthetic spectrum (reduced to the resolu- 
tion of the observed data) onto the grid of observed wavelength 
points and then calculating 

^ / ymodcl yobs \ ^ 

1 = Yl m (^ jmodcl + jobs j 

with Wi = OMfi+l + fi hB )i X i+\ - x ° hs ) where / m ° dcl is 
the (mapped) flux of the model spectrum, / obs is the observed 
flux, and A obs the observed wavelength. For each model, this 
procedure was repeated for 0.0 < E[B — V) < 2.5 in steps 
of 0. 1 to independently determine the reddenin g. For this pro - 
cedure, we used the reddening model of Carde fliet alJ 1^989). 
We then selected the models that resulted in the 3-10 lowest 
q values as the most probable parameter range for each indi- 
vidual star. The "best" value was chosen by visual inspection, 
at this point additionally considering the optical spectra to en- 
sure a consistent fit. This procedure allows a rough estimate 
of the uncertainty in the stellar parameters. Note that it does 
not eliminate systematic errors in the stellar parameters due 
to missing, incorrect or incomplete opacity sources. The com- 
parison was done for a total of 377 model atmospheres with 
solar abundances in the range 2000 K < T c s < 5000 K and 
5.5 < logg < 0.0. Together with the search range in extinc- 
tion this leads to 7539 combinations that were considered in 
the procedure. With the exception of allowing the extincti on to 
vary th is is the same procedure that was used in lLeggett et alJ 
feOOll) and lLeggett etaD ( l2002l) . 

The best fitting model has an effective temperature T e gp = 
4 200 K and a reddening of E(B - V) = 0.9. The low reso- 
lution of the data and the relative insensitivity of the spectral 
energy distribution to gravity prevent us from determining a 
value of log g, it is clear, however, that the object is a dwarf 
rather than a giant. The formal error in effective temperature is 
about ±200 K and about ±0.2 for the extinction. The low res- 
olution data also prevent detailed metallicity determinations, 
and so far only solar metallicities were considered. Overall, the 



spectral analysis results suggest a spectral type of about K7 V 
(± 2 subclasses). 

The resulting fit is shown in Fig. [3] We have applied the 
reddening to the synthetic spectrum (dotted line) in order to fa- 
cilitate the comparison without modifying the data themselves. 
All available spectral and colour information is included in the 
figure. The fit is in general acceptable, unfortunately data are 
missing in spectral regions where they would be extremely use- 
ful to test the resulting model parameters. 

A consistency check of our solution can be performed 
by comparing our measured reddening with the model of the 
Galact ic interstellar extinction constructed by lArenou et alJ 
( 1992). First we estimate the distance fro m the spectral type - 
absolute magnitude calibration of Schmidt-Kalerl ll 1 9821) . From 
their Table 13 we get an absolute magnitude of My = 8.1 for 
spectral type K7 V. With E(B - V) = 0.9 as derived above, the 
dereddened V magnitude is 15.3 (adopting R = 3.1). Thus the 
resulting distance module is 7.2, corresponding to a distance of 
280 pc. The reddening predicted from the Arenou et al. model 
and the position of 2MASSJ05 16288+260738 (I = 178.8, 
b = -6.9) amounts to E(B - V) = 0.48 ± 0.24. The scatter 
results mostly from the patchiness of the interstellar medium in 
this region. Although this value is somewhat smaller than our 
measured reddening both values agree within the error limits. 
Note that the model of Galactic extinction provides an upper 
limit of E(B - V) = 0.72 ± 0.36 for the reddening at the posi- 
tion of 2MASS J05 16288+260738. This limit results from the 
fact that stars exceeding a certain distance are above the absorb- 
ing dust layers of the Galaxy. This allows us to rule out highly 
reddened early type stars (cf. also the independent discussion 
of this aspect in Sect. l7.2l which leads to the same result). 

6. Light curve solution 

From the overall photometric data set, a subset of 23 contribu- 
tions was chosen to create the profile used for the light curve 
solution. The subsets, marked "used for profile" in Table [l\ 
were selected according to their length, the coverage in phase 
they contributed to, the filter they were taken in (=none), and 
their reliability and quality with respect to trends. In contrast 
to the data shown in Fig.|2l each of the selected data sets was 
then cleaned from suspicious points and normalised at its max- 
imum. A folded profile with 200 points, with a phase bin width 
of 0.005 units, and with phase zero set at the minimum of pri- 
mary eclipse, was then obtained from these data, and a few 
remaining clearly unreliable points were removed. 

This light curve, formed of 187 normal points (in in- 
tensity units), normalized to unity outside eclipse, was sub- 
jected toajimn^rical solution by the application of the MORO 
code JPrechsel et alJll995h . The code is based on the Wilson- 
Devinney (1971) logistical approach, but incorporates a modi- 
fied Roche model to account for radiative interaction between 
the components and uses the SIMPLEX method as parameter 
optimization algorithm. 

The solution mode was chosen such that no a priori restric- 
tion of the system configuration was imposed (equivalent to the 
original Wilson-Devinney mode 2). The total number of light 
curve parameters for a single passband curve amounts to 17. 
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Fig. 3. Measured flux calibrated optical and IR spectra (solid line), USNO-B and 2MASS colours converted to Fa (diamond 
symbols), in comparison with the best fit model spectrum (dotted). For details see text. 



Since the observed eclipse minimum depth is only moderate 
(« 16% of maximum light), as no signature of the secondary 
except its light blocking effect is evident, and because no colour 
information follows from the white light curve, solutions tend 
to be underdetermined, especially if the adjustable parameter 
set is too large. Hence it was important to use any available sec- 
ondary information from spectroscopy or stellar atmospheres' 
theory to reduce the number of free light curve parameters and 
keep some of them at fixed values. 

No information at all is available for a possible eccen- 
tricity of the orbit, since the position of the unobserved sec- 
ondary eclipse cannot be determined, and radial velocity mea- 
surements do not exist so far. Therefore circular orbits (e = 0) 
and synchronously rotating components were assumed - as 
is mostly the case in close binary systems due to their very 
short synchronization time scales. According to the late spec- 
tral type, bolometric albedos A\ and A2 were fixed at their 
usually expected values of 0.5 for convective outer layers, and 
gravity darkening exponents gi and were set to 0.32 as 
predicted by Lucy's law (1967). Linear limb darkening coef- 
ficients are poorly known for very late spectral types. F rom 
an extrapolation of the grid s of IWade & Rucinskl |l985) and 
iDiaz-Cordoves etall {H)95) at their cool ends one obtains ap- 
proximate values of x\ = 0.5 — 0.6, which were used in the 
solutions. Values of X2 (and 52) are irrelevant due to the ab- 
sence of any measurable secondary light. 

The primary effective temperature was always fixed at 
T\ = 3 000 K, typical for spectral type M5 V, since the result 
of T c ff = 4 200 K from the spectral analysis has only become 



available recently, following the February 2003 OMEGA-Cass 
observations. This choice is however not critical, because the 
light curve solution only allows to derive the temperature ratio 
T1/T2, and from a single unfiltered curve no colour informa- 
tion can be extracted. The remaining set of adjustable parame- 
ters therefore comprised inclination i, mass ratio q = M2/M1, 
secondary temperature T2, surface potentials £!i and SI2, pri- 
mary luminosity L\, and third light I3. L2 was not indepen- 
dently adjusted, but recomputed from T2 and the secondary 
surface area over the Planck law. Trial runs showed that the 
percentage of third light I3 attributable to a possible unre- 
solved field star tended toward zero (except for solution g, 
see Table [5}, so that this parameter was subsequently fixed at 
I3 = in the iterations of all other solutions. 

Convergent solutions were achieved after numerous trial 
runs with a variety of start parameter sets (start simplices) and 
different parameter increments as starting points of the auto- 
matic iteration process, which covered essentially the whole 
range of physically reasonable parameter values. For reasons 
discussed earlier the numerical process could not be expected 
to yield a single best and unique solution. Instead, for a cou- 
ple of comparably good solutions, there was no obvious way to 
qualify one of these as definitely best representation, as judged 
from the final standard deviations of normal points from the 
synthetic curves. To give an impression of the typical scatter 
of final parameters we present a subsample of 8 different so- 
lutions with the relatively best sigma standard deviations in 
Table[5] These are sorted in a sequence of increasing q values. It 
is obvious that one can identify two groups of solutions accord- 
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Fig. 4. Meridional intersection of surface and inner criti- 
cal Roche equipotentials corresponding to the nearly semi- 
detached system configuration of solution c (see Table |5Ji; the 
substellar secondary component is close to contact with its 
Roche lobe. 

ing to the value of the mass ratio: solutions a-d cluster around 
q ~ 0.10 ± 0.02, while cases e-h yielded q ~ 0.18 ± 0.01. 

A common feature of all solutions are consistent values 
of inclination (i ~ 72° - 75°), temperature ratio (T 1 /T 2 ~ 
1.6 — 1.8), ratio of radii (r 1 /r 2 around 0.9), extremely low sec- 
ondary luminosity (L2/L1 « 1 — 2 • 10~ 3 ), and similar sys- 
tem configuration. As shown for solution c in Fig.|4] which can 
be considered as representative for the group of solutions with 
q ps 0.1, the secondary is of about the same size as the primary, 
and nearly fills its Roche lobe in a close to semi-detached con- 
figuration. The photometric determination of T2 and hence the 
temperature ratio must be considered very uncertain, because 
of the missing secondary eclipse and the extreme luminosity 
ratio. 

The overall representation of the observations by the theo- 
retical light curve is very good. Figure|5](top) displays the nor- 
mal points in comparison with the synthetic curve (solid line). 
Especially the eclipse minimum is matched in detail. The stan- 
dard deviation amounts to only 7.5 mmag, which corresponds 
to the typical scatter of measurements binned to normal points. 
As shown in the bottom part of Fig. [5] most observations lie in 
a la band, and all within a 3a belt, with no apparent system- 
atic deviations. Figure[6] gives a 3-dimensional impression of 
the aspects of the system at different phase angles as viewed 
under an inclination of 74°; the configuration corresponds to 
the parameters of solution c. 

7. Discussion of alternate configurations 

The system is located only 6.9° above the galactic plane (in 
an outward direction). This implies that the reddening through 
interstellar extinction is potentially very high. Although from 
spectral observations in conjunction with detailed Galactic ex- 
tinction models many stellar spectral and luminosity types 
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Fig. 5. Top part shows the observed light curve in white light 
(dots are normal points formed by binning individual observa- 
tions to phase intervals of width 0.005) together with the theo- 
retical curve (solid line) corresponding to solution c of Table|5j 
maximum light (intensity) was normalized to unity, and phases 
were computed according to the ephemeris of Sect. 13.21 bottom 
part shows residuals of observations (in intensity units) with la 
and 3a belts. 

other than late main sequence stars can be excluded, it is also 
instructive to make use of the independent information from the 
light curve solution alone. Through the geometry of the system, 
and fundamental stellar parameters that cannot be substantially 
altered even when a star resides in a close binary system, most 
of the following alternative combinations can be excluded. This 
in turn justifies the restriction of the discussion in Sect.|5]to a 
late dwarf system. 

7. 1. Reddened giant stars 

Up to spectral types G5 or earlier, luminosity class III red giants 
have radii larger than the orbital separation dictated by the mea- 
sured orbital period and a total mass sum of the system of up to 
twice their own mass. This is illustrated in Fig.Q The solid line 
represents the orbital separation of the system components as 
a function of the total system mass for the given orbital period 
of 1^29. It is therefore a strict upper limit to the radius of any 
single component of the system. The type III giant star spectral 
types are printed at the position of their radii, once over their 
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Fig. 6. Aspects of the system at three different phases; viewing 
angle is 74°, and system parameters correspond to solution c of 
Table|U 



corresponding stellar mass and once at twice that value. These 
overplotted radii f or type III giants were taken from the mass- 
radius relation by ICoxll200(l Using the stellar mass and the 
double of it means that in between those two values all possible 
mass combinations are covered, since the roles of primary and 
secondary would simply become reversed if an even greater 
fraction of the total mass were attributed to the presumed sec- 
ondary. The first case represents the limit in which the mass 
of the companion is negligible, so that the total mass is solely 
made up of the giant's contribution, while in the second case 
the mass of the giant amounts to half of the total mass in the 
system. 

In a naive consideration, the early type giants could fit 
within the orbital separation, even if it is clear that most of the 
time they would reach well over half of the total distance. But 
although these earlier types could just about fit into the system, 
it can easily be shown that their deformation within the Roche 
potential would in all cases result in large ellipsoidal light vari- 
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Fig. 7. Orbital separation a as a function of the mass sum Mi + 
M2 of the system at the 1.29 d period. Overplotted are radii for 
type III giants; explanation see text (Sect. l7.lt . 

ations, which are not observed in the actual light curve. In 
addition, the width of the observed eclipse minimum would 
cover a much broader phase range. For these estimates, the bi- 
nary eclipse simulation program nightfall (R. Wichmann, 
Landessternware Heidelberg, Germany) which calculates syn- 
thetic light curves taking into account the distortion of the stars 
in Roche geometry was used. 

Since these considerations equally apply to luminosity 
classes II and I, and even in a much stricter form there, the lu- 
minosity class for the more luminous object in the system must 
be V or higher. 

7.2. Reddened earlier main sequence stars 

Giant stars do not fit within the prescribed orbit; but what about 
bright early main sequence stars that appear reddened by strong 
interstellar absorption? Very early main sequence stars have 
masses and especially radii similar to those of type III giants, 
so that, as above, geometry arguments can be brought forward 
to rule out a combination of two very early-type components. 
This is important, as results from stellar structure will be uti- 
lized to find physically meaningful pairs in what follows. 

The light curve solutions constrain the mass ratio, the ra- 
dius ratio and the ratio of the effective temperatures almost re- 
gardless of the absolute value of Ti. Using tabulated values for 
the masses, radii and effective temperatures of stellar and sub- 
stellar objects, the possible components making up the system 
can be constrained by requiring that both of them have param- 
eters reasonably close to those of isolated main sequence stars 
or substellar objects. The s t ellar p arameters used in the follow- 
ing were taken from ICoxl d2000h . those for su bstellar objects 
(for ag es r anging from 1 M yrs t o 10 G yrs) from Cha brier et alJ 
(hoOOh and Baraffe etalJ J2002I) . 

The ratios of effective temperatures and masses of the bi- 
nary components for the light curve solutions c and g from 
Table [5] which can be considered representative for the two 
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bulges of solutions clustering around q = 0.10 and 0.18, are 
used to find the corresponding effective temperatures and 
masses of the secondary as a function of primary mass. For 
all main sequence stars, their effective temperatures and those 
required by the two representative solutions for the secondary 
are plotted in Fig.[8] The zero-age main sequence objects and 
the youngest substellar objects (lMyrs old) are marked by 
plus signs and are connected by a thick solid line; higher age 
substellar models are also represented by plus signs, which re- 
main, however, isolated for clarity. For both solutions a row of 
squares connected by a solid line is shown. The squares corre- 
spond to the locations of the secondary in the T c g — M dia- 
gram, which follow from the temperature and mass of a main 
sequence primary using the temperature and mass ratios of the 
respective photometric solutions c and g. For each square plot- 
ted, the corresponding error estimates from the typical disper- 
sion within each of the two groups are indicated by small dots 
which represent the end points of the associated error bars (not 
drawn as full lines to preserve more clarity in this complex rep- 
resentation). For the primary a variety of spectral types between 
05 V and M8 V were considered to cover the full zero-age 
main sequence (plus signs). When inspecting this figure and 
the following graph, note that the plot scale is logarithmic so 
that offsets between curves can be much larger in regions of 
the plot that correspond to the upper main sequence than they 
might intuitively seem. 

The curve for the secondary corresponding to solution c 
only approaches and intersects the main sequence at its lower 
end and therefore excludes highly reddened hotter main se- 
quence stars as a possible primary, since the corresponding sec- 
ondaries cannot exist. The other line corresponding to solution 
g starts off close to the main sequence and comes back to it ear- 
lier than the other one. As stated above, the combination of two 
upper main sequence stars as a possible solution can be ruled 
out, because such extended stars could only reside within the 
given orbit if an appreciable distortion of the primary is allowed 
for, which would inevitably result in an easily observable ellip- 
soidal light variation. Apart from this special case for g on the 
upper main sequence, solutions (discretised in, on average, 5- 
subclass steps ! ) were elected possible whenever the error range 
for such a discrete secondary location intersected the stellar 
or substellar regime. Errors in parameter ratios are regularly 
smaller than the discretisation in spectral classes used, so the 
limits given can be regarded to be accurate to within roughly 
two subclasses. 

On the lower main sequence, the earliest possible spectral 
types for the primary in the two cases are as listed in the first 
line of Table [6] In case g, K5 and MO primaries must be ex- 
cluded. When the mass-radius relation is taken into account in 
addition to the just invoked mass-temperature correlation, these 
upper limits can be even further constrained, as will be shown 
next. 

The ratios of radius and mass for the binary components 
(also taken from Table Q together with tabulated mass-radius 
relations from the same sources as above can be subjected to 
the same procedure. The designations in Fig.|9]are analogous 
to those in Fig.[8] The results are also compiled in Table[6] The 
primary can be constrained to be of spectral type G5 or later 



for solution c, or of spectral type MO for solution g. The sec- 
ond line in Table[6]lists these limits without any additional age 
constraints that migth be present (see this discussion later). 

Combining the constraints from both Fig. [8] and Fig. [9] 
yields an overall upper limit for each solution as listed in 
Table [6] line 3. Solution c allows for a primary no earlier than 
G5, while solution g restricts possible primaries to spectral 
types no earlier than M2. Spectroscopic results strongly favour 
the group a-d solution, since the overall constraint of G5 for 
the primary spectral class is entirely consistent with the con- 
clusion in Sect.[3] For an upper limit of M2, on the other hand, 
it would be hard to claim consistency with the spectroscopy 
results. Table [6] nevertheless explores the mass range for the 
secondary in different scenarios (entries in lines 5 and 6). 

An additional constraint not taken into account so far is the 
age of the system, which for the more likely solutions a-d is re- 
stricted to below 0.01 Gyrby the mass-radius relation. A young 
system is also allowed for by the mass-temperature relation. 
However, this corresponds to a lifetime of the system in which 
the K star will not have had enough time to attain the zero- 
age main sequence yet, and will hence not necessarily have the 
ZAMS parameters assumed to deduce these constraints in the 
first place. This might well limit the overall usefulness of this 
discussion, and is a point that will have to be re-addressed once 
better data has become available for this object. 

7.3. Nearly identical components 

A serious objection to the interpretation presented so far 
emerges if the orbital period is really twice as long as assumed 
(see a similar initial ambiguity for CM Dra where this ind eed 
turne d out to be the correct interpretation in Eg gen & San dage 
Il967t) . The eclipse ephemeris given in Sect. l3.2l would then not 
correspond to the orbital ephemeris of the system any more, as 
assumed throughout the light curve analysis in Sect.[6] Hence, 
results obtained there are not applicable to the current discus- 
sion, where identical components could then produce undistin- 
guishable primary and secondary eclipses. In many binary sys- 
tems, the mass ratio is close to one, so this is not an altogether 
implausible, but from a statistical point of view highly unlikely 
configuration. This possible complication can currently not be 
resolved, since the scenario could only be conclusively ruled 
out, or corrobated, with radial velocity measurements. 

7.4. Two old white dwarfs 

A further scenario that requires consideration is a system con- 
sisting of two old and therefore very red white dwarfs (the 
effect of interstellar reddening cannot contribute significantly 
here due to the low intrinsic luminosity of these objects). It is 
however not very probable that the mass ratio in a double de- 
generate system is as low as q « 0.1. For the given period, the 
duration of eclipse would be of order 10~ 3 phase units, com- 
pletely incompatible with observations. 
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8. Conclusion 

Despite the remaining uncertainties, from the data presented 
in this paper it is plausible that the newly discovered eclipsing 
binary system 2MASS J05 16288+260738 consists of a late K- 
type (pre-)main sequence star as a primary and a substellar ob- 
ject as a secondary. For the spectral class upper limits derived 
in Sect. 17. 21 Table [6] also lists the secondary masses according 
to the mass ratios given in Table [5] for each solution. All of 
these mass values, which are close to or below the substellar 
limit of 0.075 Mq required for stable hydrogen burning, repre- 
sent upper limits. Taking into account the additional informa- 
tion available from spectral analysis which favours a spectral 
type around K7 results in a value of «0.06Mq for the sec- 
ondary's mass. In this case, only solutions a-d are considered as 
likely since a spectral type of K7 would not be consistent with 
solutions e-h. A substellar nature of the companion is there- 
fore quite likely: The unusually low mass ratios in all solutions 
make the secondary a good Brown dwarf candidate. 

This interpretation should now be checked by trying to con- 
firm the spectral classification via the detection of spectral lines 
in new high resolution, high signal-to-noise optical and/or in- 
frared spectra. These lines could then also be used to obtain ra- 
dial velocity measurements for the system which should even- 
tually provide absolute masses. 

Together with the extensive light curve available, the sys- 
tem has the potential to provide a new high-quality point for 
the mass-radius relation of the lower main sequence (or for pre- 
main sequence evolutionary tracks), and the first one obtained 
from eclipse measurements for a sub-stellar object. 
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Table 5. MORO solutions of the light curve of the eclipsing system 2MASS J05 16288+260738 
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Parameter 


a 


b 


c 


d 


e 


f 


g 


h 


q (= M2/M1) 


0.085 


0.094 


0.110 


0.122 


0.156 


0.181 


0.182 


0.184 


/. 


75? 5 


74? 3 


74? 


73? 9 


73? 5 


72? 4 


72? 9 


72? 3 


T1/T2 


1.67 


1.72 


1.74 


1.73 


1.77 


1.63 


1.83 


1.88 


n/rg 


1.11 


0.89 


0.85 


0.87 


0.85 


0.67 


0.73 


0.64 




5.138 


5.350 


5.510 


5.520 


5.598 


6.123 


6.001 


6.220 




1.961 


1.943 


1.997 


2.055 


2.182 


2.197 


2.236 


2.198 


L\ 


0.998 


0.998 


0.999 


0.998 


0.999 


0.993 


0.999 


0.999 


x\ 


0.60* 


0.60* 


0.50* 


0.50* 


0.60* 


0.60* 


0.488 


0.60* 


X2 


0.50* 


0.50* 


0.50* 


0.50* 


0.50* 


0.50* 


0.527 


0.50* 


id 


0%* 


0%* 


0%* 


0%* 


0%* 


0%* 


1.5% 


0%* 


la deviation 


0.00750 


0.00750 


0.00749 


0.00749 


0.00749 


0.00755 


0.00750 


0.00749 



"Ratio of mean Roche radii; 

^Relative luminosity Li/(L\ + L2); L2 is not independently adjusted, but recompute d from T2 and T2 ', 
c Linea r limb darkening coefficient; theoretical value for V band taken from iDiaz-Cordoves et"al] 
11995); 

d Fraction of third light at maximum; 
* fixed. 



Table 6. Interpretation of Figs.|8land|9j for explanations see Sects. 17. 2l andl8l 



possible primary spectral types for solution 


group a-d 


group e-h 


from mass-temperature relation (Fig.|8j 
from mass-radius relation (Fig.|9j 


upper limit G5 
upper limit G5 


region GO - K0, or upper limit M2 
upper limit M0 


combined constraints 

consistency with spectroscopy (K7±2) 


upper limit G5 
fully consistent 


upper limit M2 

marginally consistent at most 


resulting secondary mass [Mq] 
resulting secondary mass [Mq] for K7 


upper limit 0. 1 1 
0.062±0.01 


upper limit 0.076 
0.105±0.01 
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Fig. 8. Effective temperature as a function of the stellar mass for zero-age main sequence stars and substellar objects (plus signs) 
and their secondaries according to Table[5](squares); further explanations see text (Sect. l7.2l . 



1 00.00 F 




Fig. 9. Mass-radius relation for zero-age main sequence stars and substellar objects (plus signs) and their secondaries according 
to Table|5](squares); further explanations see text (Sect. l7.2> . 
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